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Abstract. The velocity distribution of neighbouring stars is 
deduced from the Hipparcos proper motions. We have used 
a classical Schwarzschild decomposition and also developed 
a dynamical model for quasi-exponential stellar discs. This 
model is a 3-D derivation of Shu's model in the framework of 
Stackel potentials with three integrals of motion. 

We determine the solar motion relative to the local standard 
ofrest(LSR)([/ Q = 9.7±0.3kms _1 , V Q = 5.2 ±1.0 km s" 1 
and Wq = 6.7 ± 0.2 km s -1 ), the density and kinematic 
radial gradients, as well as the local slope of the velocity 
curve. We find out that the scale density length of the Galaxy 
is 1.8 ± 0.2 kpc. We measure a large kinematic scale length 
for blue (young) stars, R ar = 17 ± 4 kpc, while for red 
stars (predominantly old) we find R ar = 9.7 ± 0.8 kpc (or 
R a 2 = 4.8 ±0.4 kpc). 

From the stellar disc dynamical model, we determine 
explicitly the link between the tangential-vertical velocity 
(ve,v z ) coupling and the local shape of the potential. Using 
a restricted sample of 3-D velocity data, we measure z , the fo- 
cus of the spheroidal coordinate system defining the best fitted 
Stackel potential. The parameter z Q is related to the tilt of the 
velocity ellipsoid and more fundamentally to the mass gradient 
in the galactic disc. This parameter is found to be 5.7± 1.4 kpc. 
This implies that the galactic potential is not extremly fiat and 
that the dark matter component is not confined in the galactic 
plane. 

Key words: Stars: kinematics - Galaxy: fundamental parame- 
ters - Galaxy: kinematics and dynamics - solar neighbourhood 
- Galaxy: structure - dark matter 



1. Introduction 



The Hipparcos proper motion measurements ( ESA, 1997 ) of 
neighbouring stars radically enlarge the size and quality of 
kinematical unbiased samples and provide a new opportunity 
to reconsider the local galactic structure and solar neighbour- 
hood kinematics. Visualization and analysis of the local veloc- 
ity field from the Hipparcos tangential velocities (Chereul et 

Send offprint requests to: O. Bienayme 
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1997 



1998) or combined with radial velocities (Figueras et 



al 

al, [19971) show the large structuring of phase space and reveal 
the process of phase mixing with time. Clearly a much better 
understanding of the local kinematics and of the galactic dy- 
namics will be obtained while the real complexity of the local 
kinematics is now apparent. 

In this paper, we re-examine the classical kinematic anal- 
ysis. We focuss this study on the two following points: 
1) it is necessary to build a dynamically coherent model of the 
kinematics since an analysis based on the Jeans equations is 
not sufficient and a more rigorous approach follows from using 
solutions of the Boltzmann equation; 2) we measure the galac- 
tic gradients (stellar density and kinematics, slope of velocity 
curve) from the local kinematic data and constrain the potential 
shape and the radial forces close to the plane. We deduce the 
flattening of the galactic potential (and mass distribution) in the 
Galaxy. This is a complementary approach to the more classi- 
cal K z problem, i.e. the analy sis of the force perpendicular to 
the galactic plane (Creze et al, 1998 ). 

In this paper, we analyse the Hipparcos proper motions, 
first using only a Schwarzschild decomposition (§2) then with 
the help of a 3D dynamical galactic model (§3) which is a 
generalization of Shu's model. Conclusions regarding the so- 
lar motion relative to the local standard of rest (LSR), density 
and kinematic gradients are drawn (§4). We finally present a 
first observational constraint on the radial potential bending at 
the solar position (§5). 

2. Kinematic data 

The Hipparcos data allow us to build a kinematic sample free 
from kinematic bias, that is the properties used to define the 
sample are not a kinematic parameter or do not explicitly de- 
pend on velocities. Survey stars, labeled S on column H68 in 



HIP catalogue (ESA, 1997), satisfy th is criteria. T hey define a 



magnitude limited sample from HIC ( |ESA, 1992[ ) magnitudes 
with dependence on galactic latitude, V < A + 1.1 sin \ b\, with 
A = 7.9 for spectral types earlier or equal to G5 and A = 7.3 
otherwise. 

A first sample has been defined with Survey stars by re- 
moving spectroscopic binaries (flag H59 set) and potential bi- 
naries in the Hipparcos solution (flag H61 set to S) and keep- 
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ing stars with relative errors on parallaxes smaller than 10% to 
avoid kinematics containing erroneous or additional motions. 
The resulting sample (11766 stars) has been divided into 19 
sets corresponding to successive ranges of B — V colour for 
main sequence stars, and a twentieth set for giants in order to 
examine the kinematic variations between the subsets. 

We have used the proper motions and parallaxes giving the 
tangential velocities to probe the local kinematics. Correction 
for gal actic rotation has been applied using Feas t & White- 
lock's ( 1997 ) values as done by Binney et al (1997). 

We also analysed a second sample in order to use the three 
components of velocity, adding radial velocity measurements 
to the tangential velocities, and expecting to obtain a stronger 
constraint on the galactic structure. However, it is known that 
radial velocities have been measured preferentiall y for high 
proper motion stars (see figure 2 from Binney et al, 1997) and, 
as a consequence, without careful selection or analysis, such 
data may introduce kinematic bias. That reason alone justifies 
the necessity to analyse firstly the Hipparcos proper motions 
without radial velocities. 

In order to analyze the full three kinematic components, 
v r , vg, v z , we have used stars from the McCormick survey 



(Vyssotsky, 1963). This sample (895 stars) is based on spec- 



troscopic surveys of identification of M dwarf s and must be 
without kinematic bias (Weiss & Upgren, 1995 ). Radial veloc- 
ities are obtained from the literature; most data were found in 



the following cat alogue s: HIC ( pSA, 1992[ ), and more recently 
i n Kam per et al ( |1997| ), Reid et al ( |1995| ), Upgren & Harlow 
( 1996). The Simbad data base has allowed for a systematic 



search and cross-identifications. Keeping only data for which 
velocity accuracy is better than 10 km s _1 and rejecting iden- 
tified spectroscopic binaries and variables, we obtain a sample 
of 309 stars. These stars have kinematics slightly hotter than 
the sample of giants. They have been combined to a sample of 
giants and red stars from HIP survey stars. 

2.1. Preliminary analysis 

To understand the main behaviour of the distribution function 
f(u, v, w) for the various stellar samples, we refer to the anal- 



ysis given by Chereul et al (1997 



1998), Gomez et al (1997) 



Figueras et al (jl997|), Binney et al fll997[ ), Dehnen ( |1998[ ), giv- 
ing both the axis ratio of velocity distributions and descriptions 
of clumps in phase space and irregularities. We have proceeded 
to a Schwarzschild distribution fitting (more exactly to a 3D 
gaussian fitting with one of the principal axes fixed to z axis) 
to the Hipparcos proper motions and (I, b) coordinates distribu- 
tionQ 



A similar analysis has been done by Dehnen & Binney (1998) 
who use a closed form to invert the Hipparcos tangential motions into 
mean solar velocity and velocity dispersions. The clear advantage of 
their method is to avoid any prior on the underlying velocity distri- 
butions (Gaussian for instance). While our classical Gaussian fitting 
give mean velocities of stellar groups close to the values they obtain, 
we obtain different results concerning the dispersions. Performing nu- 
merical simulations we find that their method introduces a bias on 



Table [T] and Figs. |l]-^-||give results for each colour bin. U is 
the velocity towards the Galactic center, V towards the Galac- 
tic rotation and W towards the north Galactic pole. Fitting with 
two ellipsoids does not improve the determination, the second 
fitted ellipsoid adjusting one of the moving group. We find the 
solar velocity relative to the local standard of rest (LSR) to be 
U Q = 9.90 ± 0.20 km s" 1 and W = 7.05 ± 0.10 km s _1 . 
The Vq component is usually taken as the limit a\ — > in 
the asymmetric drift relation (Fig. ^|). However the slope is not 
constant and the limit may range between 4 and 8 km s _1 de- 
pending heavily on the range of a\ variances selected. The 
fact that the asymmetric drift relation is not linear cannot be 
related to the a w ja u variation that is too small. The simplest 
explanations should be that the density or the kinematic scale 
lengths change gradually with populations or that the bluest 
populations are not kinematically relaxed. Such featu res are 
also marginally visible in Fig. 1 given by Mayor ( 1974 ). 

The vertex deviation is maximum for blue stars (27°) and 
for the reddest stars it fluctuates around 6°. Since red stars are 
a mixture of old and young stars with a large vertex deviation, 
it is compatible with a null vertex deviation for the old stars. 
Accurate ages help in a nswer ing this question and show a null 
deviation (Gomez et al, 1997). 



The velocity ellipsoid axis ratio a w /cr u varies rapidly with 
colour bins, increasing from 0.45 for blue stars to 0. 55 to red 
ones, similar to the feature obtained by Gomez et al ( 1997 ) by 
splitting stars by age intervals. 

2.2. Model fitting, parameter estimates and errors 

The Schwarzschild distribution in this section and the dynam- 
ical models in next ones represent velocity distribution func- 
tions f(v r ,vg,v z ) for stellar disc populations. The predicted 
tangential velocities distribution f(l, 6; Ti, Tb) as a function of 
galactic coordinates I, b and tangential velocities TJ, Xj, are eas- 
ily deduced. Stars with extreme velocities larger than 3 times 
the mean dispersions have been rejected (11685 stars remain- 
ing). This allows for the removal of possible halo stars or ob- 
jects that could not be described by our disc star modeling. All 
models are fitted to data by adjusting the model parameters 
by a maximum likelihood. This means that the fitting proce- 
dure gives unbiased estimates of the model parameters (and is 
asymptotically convergent) if the observations are drawn from 
a parent distribution corresponding to the assumed model. Er- 
rors are obtained from the second derivatives of the likelihood 
at the position of maximum. The inverse of this matrix is the 
covariance matrix of the derived parameters giving errors and 
correlations. The estimated errors have also been checked by 
repeating fits to random subsamples. 

Thus errors on parameter estimates are dominated by the 
data sampling. A much smaller source of errors comes from the 
accuracy of proper motions and distances measured by Hippar- 



estimated dispersions. The bias increases with the ratio of the mean 
velocity to the dispersion and cannot be neglected for the a v disper- 
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Table 1. Velocity dispersions (km s _1 ) and vertex deviation 
(degree) for the twenty stellar groups. Colour bins 1 to 18 have 
500 stars, the 19 th 411 and giants groups, the 20 th , 2208 stars. 
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cos. A crude estimate of such errors is obtained by considering 
the median error on tangential velocities (from proper motions, 
parallaxes and their respective errors within our sample) and 
we find 1.3 km s _1 for one velocity component. With a sample 
of 1 1 600 stars, the resulting error, for instance, on each solar 
velocity component is only 0.01 km s -1 . On the other hand, 
systematic effects, such as correlations between the Hipparcos 
solutions for close stars, exist while at large scale they must be 
null at the exception of the error related to the link to an iner- 
tial reference system. The deviation from inert ial, non-rotating 
within ±0.25 mas/yr (Kovalevsky et al, 1997 ), has negligible 
effects on the estimate of moments. 

Systematic errors or bias due to the inadequacy of the 
model are naturally the most difficult to evaluate. Only a more 
general model would help to quantify that question. For in- 
stance the non-zero vertex deviation is not considered by Eq. 1 
the axisymmetric model (a non-axisymmetric model with three 
integrals could have been considered, but in the case of a 
Stackel potential, the vertex deviation would have been the 
same for all stellar groups). Certainly the vertex deviation in- 
troduces a bias on the estimate of velocity moments if the mo- 
ments are estimated assuming a null-deviation. This will result 
in a bias on the other estimated parameters, scale length, slope 
of the rotation curve and is discussed further. 

A non-parametric approach could seem more attractive, 
avoiding the bias of parametric models and allowing for a pos- 
sibly unbiased determination of the various velocity moments. 



Fig. 1. Mean velocity ( km s x ) and vertex deviation depen- 
dences on stellar groups in colour bins. 
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Fig. 2. Velocity dispersions and ratios. 



However, in this paper, we want to avoid a direct determination 
of the velocity moments and then the use of relations drawn 
from Jeans equations since we have shown that such e quations 
lead to bias (see §3.2 and B ienayme & Sechaud 1997 , see also 
Kuijken & Tremaine, 1991 ) and that a more rigorous approach 
is acheived using solutions of the Boltzmann equation. 
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Fig. 3. Tangential velocity drift versus velocity variance. 



3. Kinematic model of galactic stellar discs 

The following analysis is based on a parametric model for 
an axisymmetric thin rotating stellar disk. Such model was 
firstly introduced by Shu (1969) and extended in Bienayme 
& Sechaud (1997) adding independent gradients for the den- 
sity and the kinematics (see also Kuijken & Dubinski, 1995). 
This model is a parametric distribution function having a two- 
integral of motion dependence, the energy E and the angular 
momentum L z . The fundamental parameters are closely linked 
to the radial density gradients of the disc stellar populations and 
also to their kinematic gradients and to the shape of the circular 
velocity curve. It has been applied to the analysis of the CNS3 



(Gliese & Jahreiss 1991 ) in order to explain the observed local 



velocity distribution (Bienayme & Sechaud, 1997). 

The main improvement of this model compared to Shu's 
model is its ability to reproduce the observed marginal veloc- 
ity distribution of stars. This is due to the introduction of three 
realistic independent parameters, the density, kinematic disper- 
sion and rotation curve gradients (see also previous investiga- 
tions of the local kinematics with dynami cally c oherent distri- 



1991, Evans & Col- 



buti on fun ctions by Kuijken & Tremaine, 
lett, |1993| and Cuddeford & Binney |1994| ). 

B ased o n an application of the Jeans equation, Fux & Mar- 
tinet ( 1994 ) deduced from the observed moments (asymmetric 
drift and velocity dispersions) that the density scale length is 
close to 2.5 — 3.0 kpc with assumptions on the closure of the hi- 
erarchy of velocity moment equations (Cuddeford & Amendt, 



1992). We must however notice that such analysis using first 
and second order moments of the distribution neglects available 
informations like the velocity distribution kurtosis and skew- 
ness. 

The galactic potential can be fitted by a Stackel potential in 
the solar neighbourhood. There the effective potential can be 
expanded in powers of (i-Ro) and z. How well can this expan- 
sion be matched by a similar one for a separable potential? For 
instance, a first order development gives a r and z separable 
potential which is then used for the estimate of the force per- 
pendicular to the plane and for the mass density determination 



in the plane. Van de Hulst (1962) and Kent & de Zeeuw (1994) 



solved the expansion problem with Stackel potentials. The ex- 
pansion includes the third order term rz 2 related to the poten- 



tial flattening. If the galactic potential were of Stackel form 
we would have an exact measure of its flattening (through the 
parameter zq, the focus of the ellipsoidal coordinate system de- 
fined in the Appendix). The galactic potential is certainly not 
exactly of a Stackel form, but probably not so far from it as 
long as like here, we consider stars within a small range of 
excursion from the solar position (3 kpc radially and 1 kpc ver- 
tically). A more general analysis (using orbit computations for 
instance) could establish accurately the coupling between the 
vertical and radial motions and its link to any potential. The 
advantage of Stackel potentials is their tractability and they are 
certainly sufficient for a first analysis. In this paper it allows us 
to establish formally the importance of the potential flattening 
on the stellar motions in the solar neighbourhood. 

3.1. 3D axisymmetric stellar disc kinematic model 

The model is a 3D extension of the 2D distribution function de- 



scribed in Bienayme & Sechaud (1997). This model is defined 
in 3D axisymmetric Stackel potentials where three integrals of 
motion are known (E, L z , I3). Being a function of these three 
integrals, f(E, L z , I3) is a stationary solution of the collision- 
less Boltzmann equation. It is also close to the Schwarzschild 
distribution for small velocity dispersions and its associated 
density is nearly exponential as are the kinematic radial dis- 
tributions. 

The expression of the distribution function is given by (see 
also the Appendix for more details): 



f(E,L z ,I 3 ) 

1 1 

^o z (L z ) 



2n{R c ) S(L Z ) 
2wk(R c ) a 2 (L z ) 

,Rc{L z y , 



exp 



E - E r 



exp 



-(- 



1) 



(1) 



and is null when L z < 0, where R c = R C (L Z ) is the radius 
of circular orbit stars with angular momentum L z and having 
the energy E circ (R c ). ft is the angular velocity, k the epicyclic 
frequency and 



E(L,) = E exp[-R C (L Z )/Rv] 



<Jr{L z ) = a , r eicp[-R c (L z )/ R ar ]. 



o z (L z ) = a , z e^p[-R c (L z )/R az ] , 



(2) 



(3) 



(4) 



in the following, we put R„^ = R ar and we define the radial 
dependence of the potential in the galactic plane assuming a 
power law rotation curve v c (r) ~ r a . 

3.2. some properties 

For a convenient range of parameters (i?s, i? CTr ), such a dis- 
tribution has effectively nearly exponential density and kine- 
matic decreases with respective scale lengths close to R^ and 
R ar (Fig.Q). Models with extremely large velocity dispersions 
cannot have small scale lenghts (i.e. large density gradients). 
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In such cases, the effective density scale lengths are different 
from the input parameter J?s and such models no longer have 
an exponential decrease over a large range of radius (Fig. Q). 
Analysing the Hipparcos data, we will not be concerned by this 
problem since the velocity dispersions involved in the data are 
small. 

The asymmetric drift (AD) relates the velocity dispersion 
of any stellar group to its tangential drift. A quite general form 



is given by Binney & Tremaine, 1987 (BT) (Eq. 4.34), valid 



in the galactic plane at z = 0, relating linearly Vad to the 
variance a 2 (for constant dispersion ratios and kinematic scale 
length that is independent of stellar types). That equation can 
be considered as a first order development with a limited range 
of validity. We have computed the asymmetric drift relation 
for stationary solutions given by our model and the results are 
plotted (Fig. ||) for different sets of parameters. The AD rela- 
tion is found to be nearly linear in the range of drift from to 
30 km s _1 . This interval is apparently the range of validity of 
the classical equation. 

Hipparcos data as well as other data used to study the asym- 
metric drift fall in this range of velocity drift. However, as seen 
in §2.1, the asymmetric drift relation is not strictly linear. Fi- 
nally for stellar populations having larger drifts, like the thick 
disc, the use of a linear relation is certainly wrong. 

The asymmetric drift is neither entirely described by the 
moment equation 4.34 of Binney & Tremaine (1987) nor by 
the more exact relations such these plotted on Fig. |5| since, in 
some range of drift, different modeling may have similar AD 
relations. Higher moments or marginal velocity distributions 
help to identify and discriminate between models. This is rec- 
ognized in Fig. U where various signatures, like drift but also 
kurtosis and skewness of the vg distributions, are identifiable. 
Differences due to a short or a long scale length are visible. In 
the case of the shortest scale length 1.8 kpc, with a r increas- 
ing, f(vg) stays symmetric and the maximum shifts regularly. 
For larger scale length 2.5 kpc, the f(vg) distribution becomes 
highly asymmetric. At higher scale length 4.5 kpc, the maxi- 
mum of the vq distributions is nearly not shifted. 

3.3. f{E,L z ,I 3 )atz = 

We give, below, the simpler expressions of the distribution 
function in the galactic plane (z = 0) (see Appendix) corre- 
sponding closer to the situation of the Hipparcos sample we 
are analysing. We have: 



f(E,L z ,I 3 ) = g(r;f,6)f ± (r,z;9,z) 
with 
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where z Q (defined in Appendix as the focus of the ellipsoidal 
coordinate system) is related to the shape of a Stackel potential. 
For a fiat rotation curve with v c (r) = v Q 
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Fig. 4. Effective density and radial velocity dispersion com- 
puted from the corresponding moments of the distribution 
function Eq. 1 (models computed with z = oo and a fiat 
rotation curve v c (r) = 1, a r . e ff are plotted in the same 
units). Four radial density distributions with R a = 2R^ and 
(To = .3 — .6 — 1.2 — 1.8 are plotted. Solar position is at 
i?o/i?s = 2.5 - 3.5 and a r (R /R^) = 0.3 for the thick disc. 
Above a critical cto value, the density scale length increases and 
is no longer close to 

This vg, v z coupling is established only in the case of a 
Stackel potential and of a quasi-exponential stellar disc. Differ- 
ent couplings would have been obtained if we had considered 
another density distribution. 

Towards z Q = oo, the potential is cylindrical and 



Towards z n 



0, the potential is spherical and 



f x = exp ^- — 



^ + 1 



2 2 ' 

a 2 'v 2 



(9) 



(7) 



Dependence on the z a parameter is easily identifiable and 
can be isolated from the other parameter dependencies. This 
is discussed in §5 and in the Appendix. This dependence links 
the velocity distribution to the potential, bringing a constraint 
on the galactic mass distribution. 

4. Model versus data: solar motion and galactic scale 
lengths 

Model to data comparison is achieved by splitting the stellar 
sample by colour intervals covering various ranges of kinemat- 
ics. The ratio between the widths of the marginal distribution in 
v r to vg is mainly linked to the slope of the rotation curve, while 
the shape of the vg marginal distribution is strongly correlated 
to the scale lengths. The exact combination of these correla- 
tions is untangled by fitting the model distribution (Eq. 1) to 
the observed velocity distributions. 

4.1. Results of model fitting to Hipparcos tangential velocities 

Results are given in Tab. ^. For each colour bin we adjust one 
distribution function (Eq. 1) with the common free parameters 
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Fig. 6. Marginal velocity distributions at solar position, for 
two models with short (full lines) and long (dashed lines) 
density scale lengths. Top to bottom with increasing veloc- 
ity dispersions. Left: radial velocity distribution for both mod- 
els, middle: tangential velocity distribution for model with 
i?s = 2.5 kpc, R a — 7.5 kpc, right: for model with i?s = 
1.8kpcR CT — 15 kpc. 



Fig. 5. Asymmetric drift relations. Tangential velocity drift ver- 
sus radial velocity variance for models with various density and 
kinematic gradients (Rq = 8.5 kpc, Vq = 220 km s _1 ). Mod- 
els with the largest density scale length present a saturation ef- 
fect on the amplitude of the tangential velocity drift. 

Table 2. Solar motion relative to the LSR and galactic charac- 
teristics deduced from the local stellar kinematics (1-cr errors). 
Best estimates of parameters are shown for two ranges of data 
(2 to 20 and 3 to 20 bins) and with models with one kinematic 
scale length R ar or with two differents scales for blue and red 
stars. Typical errors are given in the last column but do not in- 
clude the correlations. Except for uq and wq, they are raised 
by a factor ~ 1.5 — 2 when correlations are considered. 





bins 2-20 
one R ar 


bins 3-20 
one R ar 


bins 2-20 

tWO R ar 


bins 3-20 

tWO R ar 




9.5 


9.5 


9.6 


9.7 ± .3 


Vq 


3.8 


3.5 


5.2 


5.2 ±.5 


Wq 


6.7 


6.7 


6.7 


6.7 ±.2 




1.5 


1.5 


1.7 


1.8 ±.1 


Ra r 


10.0 


10.0 


9.3 & 19 


9.7 ± Ah 17 ±2 


a 


-.22 


-.24 


-.22 


-.24 ± .02 


Zo 


19 


17 


23 


24 ±10 



for all colour bins being the solar motion, Ra r , a (logarith- 
mic slope of the rotation curve) and z . Parameters not given in 
Tab. ^ are the radial and vertical velocity dispersions associated 
to each colour bin. 



We find that the determination of the tangential solar veloc- 
ity vq is sensitive to the first colour bin (youngest stars) that we 
have rejected. The density scale length is also correlated to the 
solar velocity vq. All the results favour a short scale length. 

Since the AD relation is not linear (Fig. |3|), it is very natural 
to consider that the Jfe, R„ r parameters are not the same for all 
stellar populations. The change of the ellipsoid axis ratio can- 
not explain the bending of the AD relation. In external galaxies 
the change with colour of the density scale length is only about 



20% (de Jong, 1996). However the youngest stars may have a 
smaller density scale length close to the ISM one, though it also 
depends on the SFR radial dependence. More certain, the kine- 
matic scale length of youngest populations must be close to the 
ISM kinematics with a nearly null gradient. The different rate 
of dynamical heating with radius should change the kinematics 
to the observed one for old populations. For these reasons, we 
have tried a modeling with different kinematic scale lengths for 
the bluest stars (bins 3-10) and for the reddest stars (bins 11- 
20). The likelihood is improved; we find out vq = 5.2 km s _1 
and the density scale length 1.8 kpc. We find for the bluest stars 
a nearly null kinematic gradient, in accordance to our expecta- 
tion, and for reddest stars R ar is 9.7 kpc (or R a 2 = 4.9 kpc) 
close to the Lewis & Freeman 's value (1989) R a 2 = 4.4 kpc 
based on 600 giants towards the centre and anticentre direc- 
tions. 

The comparison of the maximum likelihood (ML) of re- 
spective models allows us to determine their relative goodness 
of fit: the likelihood ratio test giving a measure of the reliability 
of adding new parameters. The likelihood-ratio statistics (Ne- 



mec & Nemec 1991L Kendall & Stuart 19671) is A = ^ where 



d is the likelihood for models based on K\ components. Mod- 
els with a larger number of components K~2 may give a better 
fit with A smaller than one. The statistical significance of A is 
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Table 3. Logarithm of likelihoods of models. Likelihood ratio tests the improvement of the goodness of fits. 



models 


In {Chins 3-20 ) 


Nb parameters 


In (Lbins 2-20 ) 


Nb parameters 


Gaussian fit 


-93993.3 


126 


-97753.3 


133 


Eq. 1 : two kinematic scale lengths 


-94090.6 


45 


-97872.4 


47 


Eq. 1 : one kinematic scale length 


-94093.2 


44 


-97887.9 


46 


gaussian fit with null vertex 


-94215.2 


108 


-97992.8 


114 



obtained by comparing —2 log A to the upper percentile of a x 2 
distribution with the number of degrees of freedom equal to the 
difference in the number of parameters between the two mod- 
els. (However, these results only hold if the asymptotic normal- 
ity and effici ency of the ML estimator are satisfied, Kendall & 
Stuart, |1967| ). Likelihoods of the discussed models are summa- 
rized in Tab. |[ Summarizing: 

1) We find that Gaussian fits are significantly improved 
when the vertex is adjusted (A ~ 220 for ~ 20 more parame- 
ters). 

2) Comparing the Gaussian fit with a non null vertex to 
Eq. 1 modelings, we see that the Gaussian gives a significantly 
improved fit (and considering they have more free adjusted pa- 
rameters). 

3) More interesting is that when Eq. 1 modelings are com- 
pared to the more similar Gaussian fits with null vertex, the 
likelihood is very significantly improved considering that the 
number of free parameters is lowered by about 65. 

4) Finally comparing the two Eq. 1 models with one or two 
kinematic scale lengths, the change of likelihoods is significant 
(as we add one supplementary parameter, the fit is improved) 
(— 21ogA=5.2 and 30 for 18 and 19 bins, corresponding re- 
spectively to an improvement with a probability of 0.98 and 
l.-5*10~ 8 ). 

Errors given in the last column of Tab. || are deduced from 
the diagonal of the covariance matrix. Existing correlations 
raise the errors. A goo d indi cator is the global correlation co- 
efficient (Eadie et al, 1971) a quantity that is a measure of 
the total amount of correlation between a variable and all the 
other variables (this coefficient being one if the variable is a lin- 
ear combination of the other variables). Most correlations and 
global correlations are small except for the following variables: 
1) a correlation and a global correlation about 0.8 between w© 
and 2) R ar has correlations of 0.3 and 0.5 with v Q and i?s 
and a global correlation of 0.8; 3) a has correlations of 0.2 and 
0.3 with vq and R ar and a global correlation of 0.5. This means 
that errors given in Tab. || must be multiplied by ~ 1.5 — 2 for 
the variables vq , R„ T , i?£ and a. 

The estimated density scale length is directly proportional 
to the adopted solar galactic radius Rq (here taken as 8.5 kpc). 
It also depends on the adopted Vq in a less direct way. We find 
that a decrease of Vq by 10 percent increases the estimate of i?s 
by about 10 percent while the slope of the rotation curve is un- 
changed, a = —0.23. Rq and Vq are not accurately known but 
they are linked by the determination of the Rq/Vq ratio (from 



the Oort's constants difference A — B), it results that i?s re- 
mains nearly unchanged by the Rq and Vb uncertainties. 

It is worth noting that we find a decreasing rotation curve 
in agreement with the ne w dete rmination of the Oort's con- 
stants (Feast & Whitelock, |l997j or Oiling & Merrifield, 1998) 



(A + B = (dV/dR) = -2.4km s^kpc^ 1 ). Considering 
the constraints given by the rotation curve (Rohlfs et al, 1988), 
a decreasing rotation curve is obtained only with small Ro and 
Vq (sm aller th an the IAU recommended values). Oiling & Mer- 
rifield ( |l998t ) find R Q = 7.1 kpc and V Q = 184 km s" 1 . The 
corresponding a, the logarithmic slope of the rotation curve, is 
then —0.10 not in excellent agreement with our local determi- 
nations. 

This difference could be partly explained by the most ev- 
ident bias of our model, i.e. the assumption of a null vertex 
deviation that leads to biased estimates of the moments of the 
ellipsoidal velocity distributions. This bias can be evaluated by 
comparing the velocity moments determined with a non-null 
vertex deviation (Tab. 1) and those determined assuming a null 
deviation. The largest axis ratio estimates are changed by 10 
to 5 percent for the six first bins and to 3 percent for the 
other bins. We quantify from Eq. 4.52 and Eq. 4.35 (BT) the 
possible order of magnitude of the bias on the slope of the rota- 
tion, a could be increased by +0. 1 to —0.14, while the estimate 
of scale lengths would be increased by 10 percent. Then the 
estimated density scale length would be 2.1 kpc. Our bias esti- 
mation cannot be applied so simply, and we will consider that 
it gives at least a more realistic estimate of errors bars, much 
larger than the Poissonian errors. 

More generally we remark that it is certainly not clear how 
correct the use of equations like Eqs 4.35 or 4.52 (BT) obtained 
by assuming an axisymmetric galaxy is when they are applied 
to observed velocity distributions that do not satisfy the sym- 
metry conditions and when the observed moments are not com- 
puted along the coordinate system. Limitation and bias in the 
application of our modeling (Eq. 1) are probably also present 
when using Jeans and the related equations. 

These results are close to previous kinematic determina- 



tions of the scale lengths by Mayor (1974) (i? s = 2.2 kpc) 



or Oblak & Mayor (1987) (R ar = 10 kpc or R a i = 5 kpc) 



and Fux & Martinet (J1994J) (i? s = 2.5 kpc) based on the local 
stellar kinematics. 

The short density scale length we have determined is also 
in agreement with the following determination from star counts 
at low galactic latitudes: with IR star counts by Porcel et al 
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Table 4. Best fit models to the 3D stellar velocity distributions. 



-Re 


Zo 


Vq 


log of 


kpc 


kpc 


km s 


likelihood 


4.0 


7.9 ±2.0 


12.4 


-34.75 


2.5 


6.9 ± 1.4 


9.5 


-31.22 


2.0 


6.1 ± 1.4 


7.8 


-30.02 


1.7 


5.7 ± 1.4 


6.1 


-29.64 


1.5 


5.4 ±1.4 


5.0 


-39.80 


1.3 


5.0 ± 1.4 


3.4 


-30.49 


1.1 


4.9 ± 1.9 


1.3 


-31.84 



( 1998 ) (2.1 kpc), Ruphy et al ( |1996| ) (2.3 kpc), Kent et al 
( |1991^ (3.0 kpc) and in visible by Robin et al ( |1992| ) (2.5 kpc). 
Higher latitude star counts (in visible) like those by Siegel et al 
(1997) may give a longer scale length. Could these results be 
biased by thick disc stars? We just mention here that the larger 
thick disc scale length (3 kpc) measured by Ojha et al (1996) 
using star counts and proper motions that allow one to separate 
accurately the stellar populations. 

Accuracy on z is extremely low since it depends on ve, v z 
coupling that is not well defined by proper motions alone. At 
2cr all z values are possible. 

4.2. 3D velocity data 

3D velocity data have been used in order to improve the deter- 
mination of z a because this parameter depends on the coupling 
of velocity components. 

For that purpose we first tried to determine z a with a subset 
from the McC survey (309 stars) but the accuracy was not suf- 
ficiently improved. We augmented this sample with HIP sur- 
vey stars with existing radial velocities. While McC sample 
is free from kinematic bias, this is not the case for Hippar- 
cos stars with published radial velocities. However the com- 
pleteness (and absence of kinem atic bias) is achieved for high 
velocity stars (Binney et al 1997 ). Since z is exclusively con- 
strained by the highest velocity stars and is unconstrained by 
stars with velocities smaller than 20 km s -1 , we expect a small 
kinematic bias. 

Combining in a single sample, McC stars (309), giants 
(1239 from bin 20) and reddest dwarfs (833 stars, B — V > 
.525) with radial velocities we find that the optimal fittings are 
achieved with a four-component model (i.e.: adding four ele- 
mentary distributions (Eq. 1) with different weights So and ra- 
dial and vertical velocity dispersions cr r! o, Cz,o)- Exploring the 
maximum likelihood space, with the parameters being free or 
setting the various parameters with values obtained in Tab. ^, 
we arrive at the following conclusions: with these data z 
and vq are correlated and may vary in the range of models 
given in Tab. [|. The accuracy and errors are determined from 
the diagonal of the covariance matrix that reflects the errors due 



to the data sampling. On the other hand, varying a in the range 
-0.3 to -0.1 (with Vq fixed) has nearly no effect on the other ad- 
justed parameters. In fact this sample alone (hot kinematic) is 
suited to constrain the parameter z but not accurately the oth- 
ers. Combining results from Tab. Qto these obtained using HIP 
tangential velocities from survey stars (Tab. ||), where the so- 
lar velocity Vq ranges between 4 and 6 km s _1 , and i?s about 
1.8 kpc, we deduce that z a is 5.7 ± 1.4 kpc. 

The zq dependence on Ro is such that zq / Rq does not de- 
pend on Rq. As a consequence the parameter A (defined in §5) 
does not depend on Rq. Concerning the Vq parameter, a change 
of 10 percent produces a negligible change of ~ 1 percent on 

Zq. 

The vertical R„ w and radial R ar kinematic scale lengths 
have been supposed equal because the data do not allow us 
to determine separately the two scale lengths. They have cer- 
tainly the same (large) order of magnitude. (The kinematic ra- 
dial scale length R ar is about 9 kpc for giants, Neese & Yoss, 



1988, Lewis et Freeman 1989, a value large compared to the 



density scale length. For younger populations the kinematic 
scale length should be even larger). 

Now considering the distribution function (Eq. 1), a de- 
crease of R a ^ (at constant R a ) has the effect of increasing the 
vertical dispersion for velocities smaller than Vb and decreasing 
for velocities larger than Vq (Figs. The apparent change on 
the isocontours would mimic a change of the z a parameter. We 
notice it also changes the velocity distribution f(vg,v z = 0) 
and then the two effects (changes on R Cw and z Q ) could be 
distinguished. Finally, modifying independently both R a and 
R am in the range 6-15 kpc, gives a range of values for fitted 
z from 3.8 to 6 kpc that is nearly about the range of the errors 
given by the maximum likelihood. 

5. Local constraints on the galactic potential 

In building the 3D kinematic model, we defined the parameter 
z Q that is a measure of the radial bending of the potential. It is 
also tightly related to 6 the vertical tilt of the velocity ellipsoid 
given by 



tan 26 = 



2al 



(10) 



Hori & Liu ( 1963 ) shown that in a Stackel potential the tilt 5 
depends only on the positions: 

2 r z 



tan 26 



(11) 



The change of the tilt close to the galactic plane (where z = 0, 
6 = and first order vertical derivatives of a r>r and ct 2jZ are 
null) is 

2 _ A(r) 



<3 z (tan2<5) = d z af. 
and for a Stackel potential (at z = 0) 

A(r) = 



r z + zi 



(12) 



(13) 
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Fig. 7. Symmetrized (axis w — —7 km s _1 ) observed f(v,w) 
marginal velocity distribution. The maximum density cor- 
responds roughly to Pleiades SuperCluster position (around 
v=— 22 km s~\ w=— 5 km s _1 ) and the innermost contour is 
certainly affected by non-stationary features that cannot be 
modeled by the stationary solution given by Eq. 1 : it is the re- 
gion where selection effects are present but it is also the region 
where data do not constrain the parameter Zq. We notice the 
egg-shape of isocontours pointing towards negative V-velocity. 

A general and approximate expression valid at z = (exact 
in the case of Stackel potentials) is given by Amendt & Cudde- 



ford ( |1991D (see also Cuddeford & Amendt, [1991 
A(r) depends on potential derivatives: 



1992) where 



A(r) 



r 2 $ r2 



3$ r + r$ r 



4r$ 2 



(r,z = 0). 



(14) 



They show that A(r) may be strongly linked to the mass gradi- 
ent in the galactic plane. For the simplest cases we have that: 

1) for a spherical potential, z a = 0, A = 1 and the velocity 
ellipsoid points towards the galactic center, 

2) while for a cylindrical potential, z Q — oo, A = and the 
ellipsoid stays everywhere parallel to the galactic plane. 

5.1. Measuring z Q 

We measure (§4.1) z Q directly from Hipparcos data using the 
distribution function Eq. 1 . Better constraints are obtained us- 
ing a sample with 3D velocities since z a is strongly tied to the 
coupling between vg and v z components (§4.2). 

Figure^ shows the f(v,w) marginal velocity distribution 
for the 3D velocity sample and Fig. ^the isocontours for three 
models: the first model is the best fit model (z — 5.7 kpc) and 



Fig. 8. Isocontours of the f(vg, v z ) distribution for three mod- 
els: best fit model (continuous lines), and best fit model chang- 
ing z to (dashed) and to oo (dotted). Isocontours are egg- 
shape, best fit and zq=0 solutions pointing towards negative 
Vg-velocity like data isocontours (Fig. Q). zo — oo solution 
(dotted) has isocontours pointing towards positive Vg -velocity. 



two similar models just changing z Q to and oo. We remark 
that these three models are strictly identical on the two axis 
f(v, w = 0) and f(v = v C i rc , w) (see Appendix), and that the 
differences are important for large velocities both in v and w. 
Outside these axis, f(v, w) depends quite exclusively on z . 
There is in fact a small dependence on a, limited by the small 
range of acceptable a values from -.3 to -.1. There is also a 
dependence on the kinematic scale lengths R ar and R a ^ that is 
null when the two scales are equal and is small when the scale 
lengths are large which is the actual situation. 

5.2. Constraints on the galactic potential 

Since the determination (z Q = 5.7± 1.4 kpc, A = 0.69± .10) is 
not yet very accurate, we have not developed detailed galactic 
mass models (uncertainties on R a and R„ m could lead at most 
to A = 0.83 corresponding to a more spherical potential). We 
mention here the main conclusions that can be drawn from the 
(z Q ; A) measure in the galactic plan e (lon ger developments are 
presented in Amendt & Cuddeford, 1991 ). 

For a galaxy with a fiat rotation curve and a full fiat mass 
distribution, A(r) = 0.5 in the galactic plane. For an exponen- 
tial disc in a fiat rotation curve potential and inside a spheroidal 
halo, A(r) ~ (the expression is valid at radii where the 
disc density remains much larger than the halo density). For a 
strictly spherical mass distribution A = 1. Combining compo- 
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nents with different fiattenings, a large range of values from 
to 1 and more are possible. 

At face value, (z Q = 5.7 kpc, A = 0.69) implies that the 
potential distribution close to the sun is not extremely flattened 
and that not all the galactic mass is confined inside the galactic 
plane. 

If the galactic halo is spherical, the radial mass distribu- 
tion in the disc exponential and the rotation curve flat, then we 
deduce from the observed A that the scale length of the mass 
distribution in the disc is about 2.7 — 3.1 kpc and is smaller 
or larger if the rotation curve is respectively decreasing or in- 
creasing. Flattening the halo is compensated by decreasing the 
mass density scale length in the disc in order to keep A to the 
observed value. However the dark matter halo cannot be flat- 
tened by more than a factor of 4 or 6, otherwise the local halo 
density is not much smaller than the disc mass density and then 
A converges quickly towards 0.5. 

More precise conclusions depend on adopted values of the 
galactic structure such as the galactic solar radius Rq and also 
the local stellar disc and ISM mass density. The solar galactic 
radius i?o may be smaller than the IAU recommended value 
(Reid, |1993| ). With R a = 7.5 kpc then A = i? /(4i? E ) is 
closer to the observed value. A smaller i?o implies that the 
rotation curve is locally decreasing (Rohlfs & Kreitschmann, 



the major axis although this is a stable configuration" (Warren 
et al, |1992|). 



1988) in concordance with our findings from the local velocity 



distribution (Tab. 2). 

When we model potentials with a flat exponential disc and 
a Miyamoto halo (with a core radius and a flattening), we do 
not find models with a strictly spherical halo and having suffi- 
ciently small A if the rotation curve is decreasing. Halos with 
a density scale height around 2 kpc at R give the observed A. 
However this result depends critically on the slope of the rota- 
tion curve and with flat or rising curve, spherical halos (density 
scale height is ~ 7 kpc at i?o) are obtained. 

We conclude that an extremely flat dark halo is rejected 
at 2cr and models favour a spherical halo or halos with a lo- 
cal scale height larger than 2 kpc at the solar galactic radius. 
There is a need for more data to bring a stronger constraint and 
to determine more precisely the galactic potential. This is also 
important because the basic measured quantities are different 
from those determined in the classical K z problem, and conse- 
quently this new method is an independent test for measuring 
the mass distribution in the galactic plane. 

We may reject with certainty some galactic mass distribu- 
tions: those corresponding to (unexpected) oblate Stackel po- 
tentials with a focus beyond the solar galactic radius and an 
ellipsoid tilt pointing in the galactic plane outwards. 

Models with A = are also rejected: they could correspond 
to cylindrical potential produced by a vertical and extremely 
prolate dark matter halo with a nearly constant mass density 
up to height of about 10 kpc. Such a halo would rotate about 
its major axis in order to have an angular momentum parallel 
to the stellar disk one. This is a confirmation of collisionless 
simulations of the dark halo that do not produce such halo con- 
figurations: "There is decided lack of halos which rotate about 



6. Conclusion 

We analyse the kinematics of nearby stars from the Hipparcos 
proper motions and determine the solar motion and the kine- 
matic and density gradients of the observed stellar populations. 
For that purpose, we have extended a dynamical model (Shu, 



1969) for axisymmetric exponential stellar discs to 3 dimen- 



sions. 

We determine the solar motion relative to the LSR to be : 
U Q = 9.7 ± 0.3 km s- 1 , V e = 5.2 ± 1.0 km s" 1 and W = 

6.7 ±0.2 km s" 1 . 

Assuming Rq = 8.5 kpc and Vq — 220 km s _1 , we de- 
termine that the stellar density scale length is short, R% = 

1.8 ±0.2 kpc. 

The nonlinear shape of the asymmetric drift is more likely 
explained assuming different R ar for young and old stars. We 
find a large kinematic scale length for blue (young) stars R ar = 
17 ± 4 kpc and for red stars (predominantly old) R ar = 9.7 ± 
0.8 kpc (or R a 2 = 4.8 ±0.4 kpc). 

These results are close to previous kinematic determina- 
tions of the scale lengths by Mayor (1974) (i?s = 2.2 kpc) or 
Oblak & M ayor (|1987| ) (R„ r = 10 kpc or R a 2 = 5 kpc), Fux 
& Ma rtinet ( 1994 ) (i? s = 2.5 kpc), and Bienayme & Sechaud 
( 1 997| ) based on the local stellar kinematics. 

We also constrain locally the rotation curve. Assuming 
v c (r) = r Q , we find a = -.24 ± 0.04. 

In the frame of Stackel potentials, we explicitly determine 
the dependence between vertical and horizontal motions. We 
show that the marginal velocity distribution of tangential and 
vertical velocities f(ve,v z ) allows us to constrain the shape of 
potential in the solar neighbourhood through the parameter z 
(focus of the coordinate system associated to the fitted Stackel 
potential) or A (Eq. 13). 

From 3D-velocity data, we measure z a = 5.7 ± 1.4 kpc or 
A = 0.69 ± 0.10. That quantity is linked to the total mass gra- 
dient in the galactic plane. It favours a spherical dark halo or 
flattened halos with a local thickness larger than 2 kpc. More 
data would allow us to refine this determination. An extremely 
flat dark halo with all the mass in the galactic plane is ruled 
out. This conclusion confirms the results obtained measuring 
K z , the force perpendicular to the galactic plane, and the lo- 
cal mass density (Creze et al, 1998), a determination based on 



vertical motions and vertical density distribution of tracer stars. 
It is consistent with preliminary results from Ibata et al ( 1998| ) 
based on the kinematics of the stream associated to the Sagit- 
tarius dwarf galaxy through the halo. 

Appendix A: disc distribution function in Stackel 
potentials 

A. 1. Stackel potentials 

A description of 3D Stackel potentials is detailed in de 



Zeeuw (1985 1. Stackel potentials are more tractable in confo- 
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cal spheroidal coordinates. The prolate spheroidal coordinates 
(A, 6, v) are related to the cylindrical coordinates (r, 0, z) by 



r2= {X + a){ V + a) md z2= (A + 7)(^ + 7) ^ 



a — 7 



7 — a 



a and 7 determine the shape of the coordinate surfaces, and 
A, v satisfies —7 < v < —a < A. Surfaces of constant A are 
spheroids and those of constant v are hyperboloids. They all 
share the same foci located on the z axis at ±z D = ±(7— a) 1 / 2 . 
In the limit z — > 00, the coordinate system is cylindrical, and 
for z Q — > it is spherical. 

Then, the potential takes the form 



h{\) - h(u) 



(A.2) 



where h is an arbitrary function (h(X) and h(v) are not defined 
on the same interval). 

Pro perties of orbits in such potentials are discussed by de 
Zeeuw ( 1985 ). They all have three integrals of motions, for in- 
stance: 



E = <5>{\, V )+ 1 -{reY + \r 2 + \z 2 



L- 



rvg 



I 3 = *(\v)-l^(r 2 + (r0) 2 ) 
27 — a 

1 . r 2 „. . 9 rzrz 

( + l)z 2 + 

7 — a 



(A.3) 



(A.4) 



(A.5) 



with 



*(A,i/) = 



2 7 — a 

(t/ + 7)/t(A)-(A + 7)/ t (i/) 
(7 - a) (A - 1/) 



(A.6) 



^3 is also the vertical energy E z when the coordinate sys- 
tem is cylindrical, i.e. z Q = 00. 

A.2. 3D disc distribution function 

In a context similar to this paper, i.e. , an analysis of a stellar 



velocity distribution, Statler (1989) used Stackel potentials and 



built a distribution function valid towards the galactic pole. He 
showed that the change of the orientation of the velocity ellip- 
soid may bias the estimated dynamical mass in the disc up to 
10 percent for samples at a distance of 1 kpc from the galactic 
plane. Analytic axisymmetric galaxy models with three inte- 
grals of motion are also described by Dejonghe & de Zeeuw 
fll988|). 



Here, we give a 3D stellar disc distribution function that 
has nearly a Schwarzschild distribution behaviour in the limit 
of small dispersions. This distribution has also a density which 
is nearly radially exponential. 



The distribution function is 

2n(R c ) E(L,) 



f{E,L z ,I 3 
1 1 



2itk(R c ) a 2 r {L z ) 
Rc{L z ) 2 



exp 



E-E r 



2ir ct z (L z ) 



exp < — ( 



o z u r J 



(A.7) 



and depends on integrals of motion. 
This can be written differently as: 

f{E,L x ,h) = f 3D (r,z;r,e,z) 

= f//(r\r,0) J- f±(r,z;9,z) 



27TC7, 



with 

f//(r;r,8) 



2fl £ 
2ttk a 2 



exp 



E 



II 



E r . 



(A.8) 



(A.9) 



where E 



II 



1/2 
2\ V r 



Vg) (when z=0, the expression of 
// / is exactly that of f 2 p, i .e. the 2D distribution function in 
Bienayme & Sechaud, 1997) and with 



fj_(r,z;6,z) = 



exp 



If! 

'2a 2 



R C [L Z 



^ 1 , i- L , 

+ i r\-,--,)h 



(A. 10) 



Properties of f 3 o are similar to /2c their density and kine- 
matics having nearly exponential radial decreases. 

The distribution function simplifies when z = 0, (i.e., 
v = —7). Adding an arbitrary constant to h leaves the po- 
tential unchanged, so without loss of generality, we may put 
/i(-7) = 0, then: 



$(A,i/ = -7) = 



h(X) 



A + 7 



*(A,i^= -7) = 



h(z - 0) 



2 zt 



and finally we obtain: 



i+(^-i) ( r2+zl 

+ [ a 2 1 \R C (L Z ) 2 + z 2 



(A. 11) 



(A.12) 



(A.13) 



A. 14) 



In the case of a cylindrical potential, radial and perpendic- 
ular motions are uncoupled and z Q — 00, so 



fi 



exp 



z 

2a 2 



(A.15) 



For a spherical potential (i.e., z a = 0) and a fiat rotation 
curve v c {r) = v Q , so R c = L z /v c and 



fx =cxp<- 



1 



1 



2 + " 2 

at at 



-)— 



(A. 16) 
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For a spherical potential (z a = 0) and a power law rotation 
curve v c (r) = v o (r/R ) a . If r — Ro, we have 



f± = exp 



at at at \v J 



2/(l+c 



(A.17) 



Finally, we remark the two following simplifications where the 
marginal distributions does not depend on z Q : 

-Firstly, for motions in and parallel to the galactic plane 
(v z — 0, z = and 7 3 = 0) we have 



f3o(r,z = 0;v r ,v g ,v z = 0) = 

f//(r;v r ,v e ) 



(A.18) 



2na z 



and secondly when z = and vg = v c (r), then R c — r and 

ho{r,z = 0;v r ,vg = v c (r),v z ) = 



f//(r;v r ,v e = v c (r)) 



1 



- exp . 
2ira 7 1 2a 2 



(A. 19) 



It is straightforward to integrate (Eqs. A 18- 19) over the ra- 
dial velocities, then similar expressions are obtained for the 
marginal velocity distributions f(vg,v z = 0) and f(vg = 
v c{r), v z ) used in §5. 

Remark: It is the requirement that the radial density (J fd 3 v) 
must be nearly exponential that makes one introduce the term 

( (j-a) — ^ ^) m ^ e seconc l exponential of Eq. A7. It is also 
this term that introduces the peculiar dependence in Eqs. A16 
or A17 instead of having the uncoupled expression of Eq. A15. 
If the density were different, the coupling in A 16- 17 would be 
different. 
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